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q , Abstract. The phase diagram of QCD at finite temperature and density and the existence of a 

critical point are currently very actively researched topics. Although tremendous progress has been 
made, in the case of two light quark flavors even the order of the phase transition at zero density is 
still under discussion. Finite-size scaling is a powerful method for the analysis of phase transitions 
in lattice QCD simulations. From the scaling behavior, critical exponents can be tested and the order 
as well as the universality class of a phase transition can be established. This requires knowledge of 
the critical exponents and the scaling behavior. We use a non-perturbative Renormalization Group 
method to obtain critical exponents and the finite-size scaling functions for the 0(4) universality 
class in three dimensions. These results are useful for a comparison to the actual scaling behavior in 
lattice QCD simulations with two flavors, as well as for an estimate of the size of the scaling region 
and the deviations from the expected scaling behavior. 
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INTRODUCTION 

In the application of lattice gauge theory to the QCD phase diagram, considerable effort 
has been spent on its successful extension to finite quark density. Nevertheless, even at 
vanishing quark density, there remains some doubt about the exact nature of the phase 
transition at finite temperature for two light quark flavors. While there is a lot of evidence 
that the phase transition is of second order yj], results by Di Giacomo and collaborators 
with staggered fermions indicate that it might be a first-order transition M>\M- 

It is difficult to establish the order of a phase transition from a simulation in a finite 
volume: The partition function remains an analytical function of the thermodynamic pa- 
rameters, no singular behavior appears and strictly speaking no phase transition occurs. 
Dynamical breaking of a continuous symmetry does not occur, either, and the analysis 
is further complicated by the necessary explicit breaking of such a symmetry. There- 
fore a careful finite-size scaling analysis is an important tool to determine the order and 
universality class of the transition. 

In QCD with two massless quark flavors, the SU(2) x SU(2) chiral flavor symmetry 
of the Lagrangian is spontaneously broken to SU(2) in the vacuum. Assuming a second- 
order transition to restore this symmetry, for massless quarks one expects this transition 
to fall into the 0(4) universality class. In the case of staggered fermions on the lattice, 
the symmetry is further reduced and one expects an 0(2) transition. This is expected if 
the QCD phase transition is dominated by the restoration of chiral symmetry. 



Results with a modified chiral lattice action (^QCD) indicate that current simulation 
volumes might still be too small and therefore outside of the finite-size scaling region [0]. 
In particular, this is relevant for the case of an 0(2) symmetry, where lattice spin model 
calculations show that the scaling region is narrow. This could account for cases where 
no second-order scaling in either the 0(2) or the 0(4) universality class is observed. 

In this context, it is very useful to analyze the scaling behavior in lattice QCD by 
comparing it to the behavior of 0(N) models. So far, mainly lattice spin models have 
been used for this purpose [0,0]. In the current contribution, we present Renormalization 
Group results for the 0(4) universality class in d = 3 dimensions. An advantage of this 
method is that we can obtain results over a very wide parameter range, which in turn 
allows for very direct comparisons to lattice results. 



RENORMALIZATION GROUP METHOD 

We employ a non-perturbative Renormalization Group (RG) method to calculate ther- 
modynamic quantities. It includes long-range fluctuations and is thus capable of describ- 
ing critical behavior. A review of functional RG methods can be found in e.g. 

The RG approach is formulated in terms of a scale-dependent effective action which 
includes quantum fluctuations between an infrared cutoff scale k and a UV cutoff A. 
The change of the potential under a change of the cutoff scale k is governed by a flow 
equation. For a specific choice of RG scheme [HSO, Q3], the RG flow equation for the 
effective potential for the O(N) model in a J-dimensional infinite volume is given by 

t ± TT ,„ m (fc¥ /2+1 1 ( (N-l) 1 \ 

5fe H ' J " (4^/ 2 r(d/2+l) \k 2 +M%{k) + &+Ml(k)J' K ) 

In finite volume, we use particular properties of our RG scheme to simplify sums over 



the discrete momenta [13]. A study of the RG scheme dependence of our results in finite 
volume is in preparation. We expand the effective potential in local n-point interactions 
around the vacuum expectation value Ob(fc): 

U k (o, n) = ao(*) + a\ (k) (o 2 + n 2 - oj(fc)) + a 2 (k) (o 2 + n 2 - o 2 (£)) 2 + ...- Ho. 

The condition 2a\(k)oo(k) =H ensures that the minimum is at (<T,7?) = (ob(£),0). 

The RG flow equation is solved numerically; input to the calculation are the values 
of the couplings at the UV scale A. For d = 3, the initial value for the minimum 
Ob (A) serves as a proxy for the temperature, (ob(A) — OQ ntlcal (A)) ~ (T — T c ). In the 
present case, we have chosen a cutoff scale of A = 1.0 GeV, which is of the order of 
the lattice cutoff (jt/a ~ 1.5 GeV) in a typical thermodynamic lattice calculation (with 
a « 0.2 — 0.3 fm) and thus appears to be a reasonable choice for a first comparison. 
Ultimately, this should be adjusted to match specific lattice results for a comparison. 



SCALING IN INFINITE VOLUME 



Critical points, such as a second-order phase transition, are characterized by a diverg- 
ing correlation length £ . The associated critical long-range fluctuations lead to universal 



TABLE 1. Results for critical exponents for 0(4) in d = 3 from RG and lattice calculations 
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FIGURE 1. Scaling behavior of the order parameter in infinite volume. Order parameter M vs. tem- 
perature t for different values of the symmetry-breaking field H (left panel), and scaled order parameter 
M/h l / s vs. the scaling variable z = t /h^ 5 l v for the same values of the field H. 



behavior where certain quantities are independent of the details of the system. Accord- 
ingly, systems can be grouped into universality classes. Close to the critical point, the 
behavior is then characterized by a small number of critical exponents specific to the uni- 
versality class. We obtain critical exponents /3, v, 8 directly from fits to the observables 
M = Oq and 1 /<§ = M a . The results are given in Tab. Q] and are in complete agreement 
with the RG fixed point analysis in ifllll . and in reasonably good agreement with spin 
model lattice results. Deviations are most likely due to the restriction to local couplings. 
Close to the critical point, the order parameter M satisfies the scaling relation 

M{t,h)=h l ' 8 f{z), z = t/h l /^ 8 \ (2) 

where z is the scaling variable and f(z) is a universal scaling function. The dimensionless 
temperature and field parameters t = (T — T c )/Tq and h — H /Hq are normalized such 
that M(t,h = 0) = (-t)P and /(0) = 1. In Fig. [Q results for the order parameter as a 
function of the temperature t for different values of H are plotted in the right panel. In 
a clear indication of scaling behavior, these curves collapse onto the universal scaling 
function f(z) when M/h 1 / 5 is plotted as a function of z (left panel). This confirms that 
the critical exponents are correctly determined and that the scaling behavior is captured. 




FIGURE 2. Finite-size scaling behavior of the order parameter at the critical temperature. Order 
parameter M vs. field h for different volume size L (left panel), and finite-size scaled order parameter 
ML?I V vs. the finite-size scaling variable hLP 5 l v for the same values of the volume (right panel). 



FINITE SIZE-SCALING 



Because universal behavior depends on a diverging correlation length % , and because 
a finite volume provides a natural infrared cutoff L, putting a system in a box is going 
to influence the scaling behavior in the vicinity of a critical point. According to the 



finite-size scaling hypothesis [14], the scaling behavior of a thermodynamic observable 
depends only on the ratio between the infinite- volume correlation length £ and the size 
of the box L. For example, for the order parameter M as a function of temperature t, the 
ratio of the order parameter Mi(t) in a finite volume of size L at temperature t and of the 
order parameter M^t) in infinite volume at the same temperature is a function of the 
ratio of the correlation length at this temperature and the volume size: 

^ = ( 3) 

Starting from this hypothesis, one can obtain universal finite-size scaling functions. 
In order to keep the ratio of the correlation length and the box size constant and thus 
keep the finite-size effect the same, one has to vary the temperature with the box size 
according to t ~ L _1 / v , since the correlation length varies as £ ~ t~ v . The situation 
is further complicated by the external symmetry breaking field h. In order to keep the 
physical behavior the same while varying the temperature, one must also keep the scaling 
variable z = t/h 1 /^) constant, and thus needs to vary h according to h ~ L~P S / V ' . 
Taking into account the infinite-volume scaling behavior, e.g. for the order parameter 
M(t,h) = h l / s f(z), one finds that the combination 

rfl v M = Q M {z,hI? 8/v ) (4) 

ought to be a universal finite- size scaling function where scaling holds. In Fig. [2[ results 
for the order parameter M as a function of the field h at the critical temperature (z = 0) are 
plotted for different volume sizes from 10 to 100 fm. In the left panel, one can clearly 
see the deviation from the infinite-volume scaling behavior. For large values of h, the 




FIGURE 3. Finite-size scaling behavior of the susceptibility at the critical temperature. Susceptibility 
X vs. field h for different volume size L (left panel), and finite-size scaled susceptibility %lJl v vs. the 
finite-size scaling variable hLP s l v for the same values of the volume (right panel). 



mass of the fluctuations is large, the correlation length is short, and deviations from the 
asymptotic infinite- volume scaling behavior occur only for very small volume size. For 
large volume size, deviations occur for small h, where the correlation length is of the 
order of the volume size. In the right panel, the finite- volume scaled order parameter is 
plotted against the scaling variable, and the curves for different volume size collapse for 
small h onto a single curve, as expected. 

In Fig. [31 results for the susceptibility % as a function of h at the critical temperature 
(z = 0) are shown. Once again, the right panel clearly shows that the finite-size scaling 
behavior is as expected, the curves for different volume size, which differ by two orders 
of magnitude, collapse onto a single curve after rescaling. The fact that scaling works so 
well also validates our results for the critical exponents, and it shows that the calculation 
method indeed includes the long-range fluctuations responsible for scaling. 

Although the results in Figs. [2] and [3] coincide quite well, the scaling behavior is obvi- 
ously not perfect. The deviations can be understood by analyzing the scaling corrections 
to the leading behavior. From an RG analysis, one can show that the corrections, which 
depend explicitly on the volume, are of the form 

L^ V M = Q^\z,hL p3/v ) + ^Q { M\zM pd/v ) + ... (5) 

The exponent co is associated with the first irrelevant RG operator at the critical fixed 
point. The coefficient functions depend only on the scaling variables z and hl?^l y . 
From the deviations, we determine CO = 0.74(4). This in good agreement with the result 



00 = 0.7338 of the fixed point analysis in the same RG scheme in Hill] . 

The main conclusion from the scaling corrections is that it might not be enough 
to analyze the leading-order scaling behavior, in particular for very small volumes. 
In the present case, where the scale is set by a UV cutoff of the order of a typical 
thermodynamic lattice cutoff, deviations are already sizable at L = 10 fm. For the 0(2) 
class, where the scaling region is narrow, a careful analysis is warranted. 



CONCLUSIONS 



In this contribution, we address the question of finite-size scaling in the O(N) univer- 
sality class with non-perturbative Renormalization Group methods. Critical exponents 
and the universal scaling functions for the 0(4) universality class in d = 3 dimensions 
in infinite volume have been obtained. The scaling behavior validates our results for the 
critical exponents and shows that our RG scheme accounts for long-range fluctuations. 

We have further demonstrated finite-size scaling behavior by implementing the RG 
scheme in a finite volume. We have calculated the universal finite-size scaling functions 
for the order parameter M and the susceptibility % for a wide range of values for the 
scaling variable z. Deviations from the leading-order scaling behavior are found to be 
consistent with the expected corrections from irrelevant operators. Due to these scaling 
corrections, an analysis taking into account only the leading scaling behavior might not 
be sufficient for small volumes used in current lattice simulations. 

Our results are suitable for a direct comparison to results from lattice QCD to check 
compatibility of the scaling behavior with the 0(4) universality class. We plan to carry 
out such a comparison, and will also extend our results to the 0(2) class. 
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